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Using cumulant expansion in Gaussian approximation, the internal quantum statistics of damped 
soliton-like pulses in Kerr media are studied numerically, considering both narrow and finite band- 
width spectral pulse components. It is shown that the sub-Poissonian statistics can be enhanced, 
under certain circumstances, by absorption, which damps out some destructive interferences. Fur- 
ther, it is shown that both the photon-number correlation and the correlation of the photon-number 
variance between different pulse components can be highly nonclassical even for an absorbing fiber. 
Optimum frequency windows are determined in order to realize strong nonclassical behavior, which 
offers novel possibilities of using solitons in optical fibers as a source of nonclassically correlated 
light beams. 



I. INTRODUCTION 

Nonclassical features of optical pulses like squeezing, 
sub-Poissonian statistics and entanglement have been of 
increasing interest for optical communication and mea- 
surement techniques |^-|^. From classical optics it is 
well known that in fibers photons can fly in soliton-like 
pulses over long distances, provided that they intense 
enough and the fiber nonlinearity can compensate for the 
dispersion- assisted pulse spreading (|J^. Yet solitons are 
not rigid but very lively nonclassical species, because of 
the quantum noise associated with the nonlinear dynam- 
ics. Since already single-mode radiation becomes non- 
classical under the influence of a Kerr nonlinearity [Q, 
solitons can be expected to exhibit not only single-mode 
nonclassical properties but also nonclassical internal cor- 
relations. 

The study of quantum solitons has been also motivated 
by a number of possible applications . Solitons may 

be used for realizing efficient and reliable sources of pulses 
with intensity fluctuations below the shot noise level. 
Nondestructive high-precision optical switching devices 
and logic gates may be implemented applying the concept 
of quantum non-demolition measurement to the collision 
interaction of solitons [p|-[rO|. Suppression of noise in soli- 
ton communication lines and usage of solitons in 
quantum information processing Q may be other poten- 
tial applications. 

Nonclassical properties of optical solitons have been 
detected in a number of experiments. In soli- 
ton squeezing is measured using homodyne detection. 
Photon-number squeezing of the spectrally filtered soli- 
tons is measured by direct detection in ||^|4j. In the 
experiment in with asymmetrical (10/90) fiber-loop 
interferometer photon-number squeezing up to 6.0dB is 
achieved. Recently the experimental studies have been 
extended to internal spectral photon-number correlations 
associated with narrow bandwidth soliton components 



Disregarding losses, the dynamics of quantum solitons 
in optical fibers may be described by the operator-valued 
nonlinear Schrodinger equation, which can be solved em- 
ploying the Bethe ansatz [0 . The method was used suc- 
cessfully to calculate field and intensity correlations in 
the space-time domain |lj,|l^. In practice however, the 
quantum nonlinear Schrodinger equation must be mod- 
ified by supplementing it with further terms in order to 
include effects such as absorption, third-order dispersion, 
and Raman scattering. Various methods of solution have 
been developed and successfully applied p0|-|3^. 

Here we follow and use cumulant expansion in 
Gaussian approximation for studying spectral photon- 
number squeezing and nonclassical photon-number cor- 
relation of different spectral components of damped soli- 
tons. We consider both narrow-bandwidth spectral pulse 
components ||l^ and pulse components with finite spec- 
tral bandwidth selected by optimally chosen square band- 
pass filters [Q. The numerically obtained results illus- 
trate the possibility to use optical solitons for genera- 
tion of squeezed light and nonclassically correlated light 
beams also in the presence of absorption. Moreover, the 
results reveal that absorption can improve, under certain 
circumstances, the nonclassical features. 

The article is organized as follows. Section || outlines 
the used basic concept. The relations necessary for study- 
ing nonclassical correlations are given in Sec. |l|, and the 
numerical results are reported in Sec. IV. Finally, a sum- 
mary and some concluding remarks are given in Sec. 0. 



II. BASIC CONCEPT 

Let us first give a brief outline of the used concept for 
describing the propagation of a quantized optical pulse 
through an absorbing fiber with second order dispersion 
and Kerr nonlinearity (for details, see and references 
therein). The spatio-temporal pulse evolution in a co- 
moving reference frame is described in terms of slowly 
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varying bosonic operators a(x, t) satisfying the commu- 
tation relation 



(1) 



[A is the effective cross-section of the fiber]. The un- 
damped motion is governed by the Hamiltonian 



H = HA dx 



(5^st) (a^a) 



(2) 



where the constant x is related to the third-order suscep- 
tibility x'"^^ as 



X 



(3) 



[er, relative permittivity at the carrier frequency Uc', 
fee, carrier wave number; Wgr, group velocity; o;^^-' = 
(Puj{k) /dk^\k=k^]- Note that solitons can be formed ei- 
ther in focusing media with anomalous dispersion (x < 0, 
tj'-^-' > 0) or in defocusing media with normal dispersion 
(X>0, a;(2)<0). 

The damped motion is treated on the basis of the mas- 
ter equation 



ih dtp = [H,p] + ijLp, 



(4) 



where 



jLp — jhA / dx 



Ntii (2a) pd — pdd^ — dd^ pj 



+ [Nth + 1) (2apat - pd'^d 
with 7 being the damping constant, and 



d'^dp) 



exp 



(5) 



(6) 



(fee - Boltzmann constant, T - temperature). The model 
applies to pulses longer than 1 ps, otherwise the influence 
of additional effects such as Raman scattering and third 
order dispersion must be taken into account. 

The master equation (Q) is converted, after spatial dis- 
cretization, into a pseudo-Fokker-Planck equation for an 
s-parametrized multi-dimensional phase-space function, 
which is solved numerically using cumulant expansion in 
Gaussian approximation. The initial condition is realized 
by a multimode displaced thermal state, without internal 
entanglement, and it is assumed that the field expecta- 
tion value corresponds to the classical fundamental soli- 
ton. In the numerical calculation, the coordinates x and 
t are scaled by the initial pulse width xq and the disper- 
sion time — x-^\uj^^'>\^^ respectively. The calculations 
are performed on a grid of 200 points with discretization 
step Ax = 0.1 xq for an initial photon number of the pulse 
of 2n = 2|w^^''^/(xa:^o)| = 2-10^ and a reservoir photon 
number of Nth = 10~^^. The value of A'th corresponds to 
a carrier wavelength of Ac = 2nc/LUc =1.5 p,m of the pulse 
in vacuum and a temperature of T = 300K [see Eq. (||)]. 
Thus the pulse is initially prepared in a displaced thermal 
state that is almost a coherent state. 



III. NONCLASSICAL SPECTRAL 
CORRELATIONS 



In order to study spectral properties, bosonic operators 
in the Fourier space are introduced, 



d{uj, t) 



1 



/2ti 



'a(x, t). 



(7) 



The operator Ni of the number of photons in a frequency 
interval i.e., Q.i<Lo <U,[, reads 



Ni^Al du d^{uj,t)d{uj,t), 



(8) 



and the photon-number variance of the beam associated 
with this frequency interval can be given by 



(ATVf) = (TVf) - (7V,)2 = (: AiV^) + (TV,), 



(9) 



where : : introduces normal ordering. When the inequal- 
ity 



(lAiV^:) <0 



(10) 



is valid - a condition that cannot be satisfied within the 
classical noise theory - then sub-Poissonian statistics are 
observed, i.e., (ANf) < (N,). 

Let us now consider two nonoverlapping frequency in- 
tervals (fli, il'j), 1 = 1,2. A measure of the mutual correla- 
tion of the photon-number variances in the two frequency 
intervals is the correlation coefficient 



m2 



(AN1AN2) 
(AiV2)(AiV| 



(11) 



(note that {AN1AN2) = (: AA^i A7V2 :) for nonoverlapping 
frequency intervals). From the Cauchy-Schwarz inequal- 
ity 

(A7Vf)(A7V|) - (ATViAiVa)^ > (12) 

it follows that the absolute value of the correlation coef- 
ficient is limited to the right. 



1^2! < 1. 



(13) 



Obviously, the correlation is nonclassical, if the inequality 
( |T^ ) is violated for normally ordered quantities, i.e.. 



{■.AN^:){:AN^:) - {■.ANiAN2:f < 0. 



(14) 



In order to give a quantitative measure of the degree of 
nonclassical correlation, we define the generalized corre- 
lation coefficient 
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{:AN!:){:ANl 



■.AN1AN2: 



{AN^){AN^)-{:ANf:){:ANi:) 



?7ll»722 



|1 - 77ll»722| 



where 



■.AN? 



Vii = 



{ANf 



< 1. 



(15) 



(16) 



Note that r]ii is negative (positive) for sub-Poissonian 
(super-Poissonian) statistics. Nonclassical correlation is 
observed if 



fju < 0, 
and it is easily proved that 

?7i2 > -1- 

From Eq. ( p^ ) it is easily seen that when 
{■.ANl.){:ANl.) <0, 



(17) 



(18) 



(19) 



then the correlation of the photon-number variances is 
nonclassical. Note that for Gaussian quantum states the 
criterion ( p^ is closely related to that used in ||l6|| . 

Similarly, the mutual correlation of the numbers of 
photons in the two frequency intervals can be considered, 



T12 



{N1N2 



{hul < 1)- Nonclassical correlation is realized if 



)-(:7Vi7V2:)2<0, 



or equivalently 



with 



T12 = 



T12 < 0, 



2 •/ 



■.N1N2 



■.Nl.){:Nl: 



(20) 



(21) 



(22) 



(23) 



(7^12 > — 1)- Note that for the photon numbers an inequal- 
ity analogous to (|l^) is not valid, since the average of the 
normally ordered square of a photon number cannot be 
negative. 

The spectral photon-number statistics can be deter- 
mined using a setup shown in Fig. 0. The scheme in 
Fig. |l|(a) is used in (and with slightly modified de- 
tection in ||l^ ) for measuring spectral photon- number 
squeezing. The scheme can also be used to reconstruct 
correlations between different spectral components from 
the measured data |l^. However with regard to direct 
correlation measurement, a scheme of the type shown in 
Fig. ^6) may be better suited for that purpose. Omitting 
the detectors, the schemes can also be used for select- 
ing from the original beam partial beams that are highly 
nonclassical. 



IV. RESULTS 



A. Narrow bandwidth intervals 



In Fig. H the temporal evolution of the coefficient 7711, 
Eq. (p6|) , is plotted for a narrow bandwidth mid-interval 
(rii,r2'i) — {—Auj,Auj), Auj/luq ^ 1, and various values 
of the damping parameter 7 (here and in the following 
frequencies are scaled by = ^^)- The strongest 
sub-Poissonian effect is observed for 7id = 0.03 at t^-^ lOtd- 
Surprisingly, the sub-Poissonian statistics of the narrow 
bandwidth mid-component can be stronger for a damped 
soliton than for the undamped one [p4[ . The temporal 
evolution of the coefficient 7711 as a function of the fre- 
quency Lu [for a narrow bandwidth interval = 
{uj — Auj, uj + Au))] is plotted in Fig. ^ for (a) 7 = and 
{b) 7td = 0.03. From the figure it is seen that the inter- 
ference pattern which is typically obtained in the limit 
of vanishing absorption [Fig. ||(a)] is not observed for an 
absorbing fiber [Fig. |3|(6)]. Note that in the limit of van- 
ishing absorption the results obtained here are in good 
agreement with those obtained by stochastic simulations 
within the frame of the positive P representation |33j . 

Disregarding absorption and solving the cumulant evo- 
lution equations given in ]3^ in a linearization approx- 
imation, the solution can be given by eigenfunction ex- 
pansion, the expansion coefficients being determined by 
the chosen initial condition. The result of superposition 
then corresponds to an interference pattern like that in 
Fig. ||(a). Obviously, the (phase-sensitive) terms that are 
superimposed do not respond uniformly to absorption, 
so that the internal coherences responsible for the in- 
terference pattern can be destroyed at least in part. In 
particular, near the center of the spectrum (lu — > 0) the 
super-Poissonian peaks are fully suppressed and in place 
of them sub-Poissonian statistics is observed. Accord- 
ingly, the super-Poissonian side-band formation that ap- 
pears with increasing frequency is more uniform for an 
absorbing fiber than for a nonabsorbing one. 

To get an insight into the correlation between spectral 
components at different frequencies tui and lu2 [for nar- 
row bandwidth intervals (fi^, Q'^) — (uji — Alu, uji + Aw), i 
= 1,2], we have plotted the correlation coefficients 7712 
[Eq. jp). Fig. I, 7712 [Eq. dH), Fig. I, and ^2 [Eq. (H), 
Fig. |6[ as functions of uji and uj2 for 7 = and 7id = 0.03, 
restricting our attention to the two characteristic propa- 
gation times t — 5td and i = 10 (cf- Fig. ^. Figures ^(a) 
and ^6) reveal that the super-Poissonian side-band com- 
ponents observed at the propagation time t = 5 id (Fig- 1^) 
arc strongly positively correlated with respect to the 
photon-number variance, whereas relatively strong nega- 
tive correlation between the super-Poissonian side-band 



3 



components and the sub-Poissonian mid-components is 
observed - effects which can be found for both a nonab- 
sorbing fiber, Fig. ||(a), and an absorbing fiber, Fig. ^(&). 
The notably reduced correlation observed at the propa- 
gation time f = lOid for a nonabsorbing fiber. Fig. ^(c), 
may be viewed as an interference effect in a similar sense 
as mentioned above. From that it can be understood that 
absorption can enhance correlation - a surprising effect 
which can be seen from Fig. ^(d). 

From Fig. [s] it is seen that nonclassical correlation 
of the photon-number variance appears between side- 
band super-Poissonian components and - in agreement 
with the condition - between the sub-Poissonian 
mid-component and the super-Poissonian sideband com- 
ponents. The relatively strong nonclassical correlation 
?7i2 ~ — 7.3- 10^"^ observed for a nonabsorbing fiber at 
the propagation time i = 5td [Fig. ^(a)] reduces to 
^12 ~ -8.6- 10-4 at t = lOid [Fig. |(c)]. The behavior 
of 7712 for an absorbing fiber may be quite different, 
as is seen from comparison of Figs. ||(a) and |5|(c) with 
Figs. and|^(d), respectively. Whereas at i = 5td the 
effect of nonclassical correlation is weaker, fji2 ^ —3-10^^ 
[Fig. ||(5)], a stronger effect is observed at t = 10 1^, 
^?i2~ -4.8-10-3 [Fig. The latter reflects the en- 

hanced correlation between the mid-component and the 
sideband components as shown in in Fig. ^(d). We are 
again left with the surprising fact that absorption can 
enhance nonclassical behavior. 

Figure 1^ shows that the photon- number correlation can 
be quite different from the correlation of the photon- 
number variance. Contrary to the photon-number vari- 
ance, there is no nonclassical correlation of the photon 
numbers of the sub-Poissonian mid-component and the 
super-Poissonian sideband components, since for the pho- 
ton number an inequality of the type (19) cannot be 
valid. Nevertheless, the minimum value fi2 ^ —0.02 ob- 
tained for a nonabsorbing fiber at the propagation time 
t — 5td [Fig. ^a)] is smaller than the minimum value 
7712 ~ -7.3- 10-3 [Fig. 11(a)], and hence a stronger non- 
classical correlation is observed for the photon number 
than the photon-number variance. Moreover, the effect of 
nonclassical correlation observed at the propagation time 
i = 10 id is weaker than that at t — bt^ for both a nonab- 
sorbing fiber [Fig. ^(c)] and an absorbing one [Fig. ||((i)] 
[contrary to the correlation of the photon-number vari- 
ance. Fig. ^(6) and {d)]. 



B. Finite bandwidth intervals 



Let us return to the coefficient 7711. Allowing for fi- 
nite frequency windows, we consider a frequency interval 
(fii, ri'j^) = (— n, f2) [for the scheme, see Fig. 0(a)] and op- 
timize n such that the sub-Poissonian components dom- 
inates the signal and rjn becomes minimal for all times 
(Fig. 0). From Fig. 0(d) it is seen that for a nonabsorbing 
fiber the absolute minimum of 7711 that is attainable in 



this way is 7711 ~ — 3.5 (i~4.5fd), which corresponds to 
a Fano factor of Fi = {ANf)/(A) = (1 - 7711)-! ~ 0.23 
or ~ 6.4 dB squeezing (cf. [^,Q). In that case ^ 80% 
of the photons of the initial, full pulse contribute to the 
filtered pulse [Fig. 0(c)]. It is further seen that the os- 
cillating behavior of 7711 as a function of the propagation 
time is damped out owing to absorption such that for cer- 
tain propagation times stronger sub-Poissonian statistics 
can be observed for an absorbing fiber than a nonab- 
sorbing one. This effect is of course a consequence of the 
behavior of th e nar row bandwidth components as ad- 
dressed in Sec. IV A| . The observed discrepancy between 
[Fig. 0(a)] and flQ [Fig. 0(5), for the definition of fip, 
see Fig. 0], e.g., for "/t^ = . . . 0.005 at i > 9 td, obviously 
results from the correlation between different frequency 
components (cf. ||l3| , ^ ). 

To obtain strong nonclassical correlation between dif- 
ferent beams, we first consider two (with respect to the 
center of the spectrum) symmetric frequency windows 
and (f72, ^^2) = (— ^^i, — f^i) [for the scheme, see 
Fig. |l|(6)] and optimize ili and fl'^ such that 7712 (Fig. 0) 
and fi2 (Fig. |9|) become minimal. From inspection of 
Fig. 0(d) the absolute minimum of 7712 is realized for a 
nonabsorbing fiber (7712 ~— 0.48 at t~4.5td)- The effect 
of absorption is again seen to damp out the oscillations 
of 7712 (as a function of the propagation time) such that 
at certain propagation times [jtd = 0.01 . . .0.02, t/t^ = 
6 ... 14 in Fig. 0(d)] stronger nonclassical correlation can 
be achieved for an absorbing fiber than a nonabsorbing 
one. On the contrary. Fig. 0(d) reveals that dissipation 
always reduces the nonclassical photon-number correla- 
tion, i.e., for any propagation time the lower bound of the 
coefficient fi2 is realized in the limit of vanishing absorp- 
tion. Note that the strongest nonclassical photon-number 
correlation (fi2 = —0.32 at t ~ 3.5 td) is weaker than the 
strongest nonclassical correlation of the photon-number 
variance. 

For weak absorption {'jtd = . . . 0.005) from Figs. 0(a) 
and (6) a (quasi-)periodic change of the frequency in- 
terval (r^i, C l't ) is seen, in agreement with the results 
in Sec. IV A, At the early stage of propagation the fre- 



quency interval is essentially determined by the nonclas- 
sically correlated super-Poissonian sideband components. 
Later [at t^St^ in Figs. 0(a) and (5)] it shifts towards 
the mid-frequency, which indicates the increasing weight 
of the sub-Poissonian mid-components. The shift back 
to the sideband components obviously results from the 
(quasi-)periodic formation of super-Poissonian (and sub- 
Poissonian) components in the center of the spectrum [cf. 
Fig. 0(a)]. The fraction of the number of photons in each 
beam relative to the number of photons in the initial, full 
beam changes from less then 5% for sideband frequency 
windows to about of 40% for near mid- frequency windows 
[Fig. 0(c)]. With increasing absorption only the first shift 
of the frequency window from the sideband to the central 
part is observed. In contrast to 7712, minimization of fi2 
always requires sideband frequency windows [Figs. 0(a) 
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and (b)], where each beam contains less than 10% of the 
initial number of pulse photons [Fig. ^(c)]. 

Whereas symmetric windows are expected to be best 
suited to realization of minimal fi2 , from Fig. |^ it is sug- 
gested that asymmetric windows may be more suited to 
realization of minimal 7712 . This is fully confirmed by the 
calculation. The symmetric windows shown in Fig. ^ in- 
deed yield the smallest value of fi2. A comparison of 
Fig. with Fig. |l^ reveals that 7^12 can be reduced 
if asymmetric windows are used. In particular, it is seen 
that the absolute minimum of 7712 observed for a nonab- 
sorbing fiber can be reduced to 7712 ~ —0.71 {t^5.5td). A 
more detailed analysis is given in Fig. |Tl|. It is seen that 
only at a early stage of propagation symmetric and asym- 
metric frequency windows yield equal values of 7712. Obvi- 
ously, asymmetric frequency windows take better account 
of the nonclassical correlation between sub-Poissonian 
mid-components and super-Poissonian sideband compo- 
nents in order to reduce 7712. The price to be paid are 
the unequal photon numbers of the two beams, since the 
number of photons in the spectral interval around the 
center is substantially larger than that in the sideband 
interval, as is seen from the figure. Hence for generation 
of (with respect to the photon-number variance) nonclas- 
sically correlated beams, the frequency windows should 
be chosen such that an optimum between nonclassical 
correlation and beam intensities is observed. 



V. SUMMARY AND CONCLUDING REMARKS 

We have studied the internal quantum statistics of fun- 
damental solitons in absorbing Kerr media, applying mul- 
timode cumulant-expansion techniques and solving the 
resulting evolution equations numerically in Gaussian ap- 
proximation. In particular, we have calculated the tem- 
poral evolution of the photon-number variance, its corre- 
lation, and the photon-number correlation for various fre- 
quency windows. The formation of super-Poissonian side- 



band components with nonclassically correlated photon 
numbers and nonclassically correlated photon-number 
variances may be regarded as being a typical signature 
of the quantum nature of a soliton pulse. 

Since a soliton pulse is a highly involved multimode 
field, interference effects can essentially determine its 
properties. It is worth noting that absorption can damp 
out interferences that are destructive with respect to non- 
classical features, such as squeezing and the nonclassical 
correlations mentioned, so that absorption surprisingly 
improves, under certain conditions, these nonclassical ef- 
fects. The calculations show that for a nonabsorbing fiber 
propagation distances of 3.5 . . . 5.5 dispersion lengths are 
best suited for detecting the nonclassical features. De- 
structive interferences are observed at distances of 8 ... 12 
dispersion lengths. At these distances the best values for 
photon-number squeezing and nonclassical correlation of 
the photon-number variance are achieved for an absorb- 
ing fiber (7td'-0.02). 

Using appropriately chosen spectral windows, a soli- 
ton pulse can serve as a source of nonclassically corre- 
lated light beams. We have considered both symmetric 
and asymmetric windows and optimized them such that 
the filtered beams are maximally nonclassically corre- 
lated with regard to the photon number and the photon- 
number variance. It should be pointed out that in prac- 
tice a number of effects such as Raman scattering and 
third-order dispersion should be included in a refined 
theoretical model. Finally, inclusion in the calculation of 
non-Gaussian effects for weak absorption has been a chal- 
lenge. 
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FIG. 1. Experimental setup for studying the spectrall y r esolved quantum statistics of soliton pulses. The scheme (a) is 
suited for measurement of the photon-number statistics whereas the scheme (fo) allows one to perform correlation 

measurements. If the detectors Di and D2 are omitted, the schemes can be used for selecting from the original beam partial 
beams that are highly nonclassical. 
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FIG. 2. The temporal evolution of the coefflcient 7711, Eq. is shown for the spectral interval (ni,r2'i) = (— Aw, Aoj), 

Auj — O.OSTTtJo, and various values of the damping parameter 7. 
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(a) 7 = 
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FIG. 3. The temporal evolution of the coefficient r;ii, Eq. (p^, is shown as a function of the frequency lo for the spectral 
interval (f2i, Sl'i) = (a; — Aw, lo + Aw), Aoj = 0.057rwo, and the damping parameter 7 = 0.03/td- Note that jyn is an even function 
of UJ. The thick solid line from the bottom to the top defines the frequency f2o [used in Fig. such that 7711 > for \u)\ > fio. 
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(a) 7 = 0, t — 5td 
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(b) 7td = 0.03, t = 5td 
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[d] 7td = 0.03, t = 10td 

1^2/ 



3 
2 
1 

-1 
-2 
-3 




3 



^12 



I 



0.28 
0.24 
0.2 
0.16 
.12 
.08 
. 04 


-0 . 04 
-0 . 08 
-0 . 12 



1 2 

FIG. 4. The correlation coefficient r;i2, Eq. (^ij), is shown as a function of the frequencies cui for the spectral intervals {Qi, H^) 
= {LJi-Auj,uJi + ALu), Alj = O.OSttijo (j = l,2), and (a) 7 = 0, f = 5td, (&) 7 = 0.03/td, t = 5td, (c) 7 = 0, f = 10td, (rf) 7 = 0.03/fd, 
t = 10fd. 
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(c) 7-0, t^lOta 
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FIG. 5. The coefRcient 7712, Eq. (|l5|), is shown as a function of the frequencies Ui for the spectral intervals {Qi,Qi) = 
ALu,uji + Alj), ALj = Omnujo (i = l,2), and (a) 7 = 0, t = 5td, (b) 7 = 0.03/td, f = 5td, (c) 7 = 0, t = 10td, (d) 7 = 0.03/fd, t = 
lOtd- Negative values of 7712 indicate nonclassical correlation of the photon-number variances at chosen frequencies ui and 0^2. 
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FIG. 6. The coefficient fi2, Eq. (^3|), is sliown as a function of the frequencies cui for the spectral intervals = 
{(jj,-ALu,uji + Auj), Alu^0.05tvujo (i = l,2), and (a) 7 = 0, t = 5td, (&) 7 = 0.03/fd, t = 5td, (c) 7 = 0, t = 10td, (d) 7 = 0.03/fd, 
t = lOtd- Values fi2 > 0.02 are not depicted. Negative values of fi2 indicate nonclassical correlation of the photon numbers at 
chosen frequencies uji and W2 . 
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FIG. 7. Optimization of the square bandpass filter [Fig.0(a)] for realization of maximum sub-Poissonian statistics. The 
temporal evolution of (a) the spectral interval (fii, il'i) = ( — fl), (c) the number (iVi) of photons ({N) is the total number of 
photons in the initial pulse), and (d) the minimum value of the coefficient 7711, Eq. (^), for that interval are shown for various 
values of the damping parameter 7. For comparison (b) , the frequency Qq defined in Fig. |^ is shown. 
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FIG. 8. Optimization of two (with respect to the center of the spectrum) symmetric square bandpass filters [Fig. ^(6)] for 
reahzation of strong nonclassical correlation of the photon-number variance between the selected beams. The temporal evolution 
of (a), (6) the spectral interval (f2i, fi'i), (c) the number (iVi) of photons in a beam {{N) is the total number of photons in the 
initial pulse), and (d) the minimum value of the coefficient 7712, Eq. (|l^), are shown for various values of the damping parameter 
7- 
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FIG. 9. Optimization of two (with respect to the center of the spectrum) symmetric square bandpass filters [Fig. 0(&)] for 
reahzation of strong nonclassical photon-number correlation between the selected beams. The temporal evolution of (a), (6) the 
spectral interval Q'l), (c) the number (iVi) of photons in a beam {{N) is the total number of photons in the initial pulse), 
and (d) the minimum value of the coefficient fi2, Eq. (p^), are shown for various values of the damping parameter 7. 
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FIG. 10. The temporal evolution of the minimum value of the coefficient 7712, Eq. (^5|), attainable for two beams obtained 
from asymmetric (nonoverlapping) frequency windows 0,'^), i — 1,2, is show for various values of the damping parameter 7. 
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FIG. 11. Optimization of two (with respect to the center of the spectrum) asymmetric square bandpass filters [Fig. ^(6)] 
for realization of strong nonclassical correlation of the photon-number variance between the selected beams. The temporal 
evolution of the spectral intervals {Qi,Qi), i = 1,2, the numbers (Ni) of photons in the beams {{N) is the total number of 
photons in the initial pulse), and the minimum value of the coefficient 7)12, Eq. (^), are shown for (a) 7 = 0, (fe) 7 = 0.01/td, 
(c) y — 0.03/td- For comparison, the dotted line shows the minimum value of 7)12 obtained for symmetric windows as given in 
Fig. ^, and the dashed line shows the minimum value of fi2 as given in Fig. |^. Note that the symmetric windows in Fig. ^ are 
best suited to minimize fi2. 
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